Itinerant metamagnetism induced by electronic nematic order 
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It is shown that metamagnetic transition in metals can occur via the formation of electronic 
nematic order. We consider a simple model where the spin-dependent Fermi surface instability 
gives rise to the formation of an electronic nematic phase upon increasing the applied Zeeman 
magnetic field. This leads to two consecutive metamagnetic transitions that separate the nematic 
phase from the low-field and high-field 'isotropic' metallic phases. Possible connection to the physics 
of the bilayer Ruthenate, Sr3Ru207, is also discussed. 
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Introduction; In terms of structural properties of 
matter, liquid crystals of molecules have intermediate 
forms between the solid and liquid phases. Common ex- 
amples of these 'soft' forms of matter include the smec- 
tic and nematic phases that have spontaneously broken 
unidirectional translational and rotational symmetries, 
and broken rotational symmetry, respectively yj. It is 
much less appreciated, however, that strongly interacting 
electron systems can give rise to electronic liquid crys- 
tal phases where the electronic structure rather than the 
structural properties of the system is 'soft' ^^j. For in- 
stance, in two-dimensional electron systems, the Wigner 
crystal and homogeneous electron liquid phases arise in 
the small and high density limits. One can ask whether 
an intermediate form of electronic structure arises in this 
case. The discovery of quantum Hall nematic phase in 
two-dimensional electron gas (2DEG) in high Landau lev- 
els clearly shows that the answer is affirmative. In this 
system, the broken rotational symmetry of the nematic 
phase is demonstrated in a series of transport measure- 
ments 0. 

When electrons reside in the lattice, broken transla- 
tional and rotational symmetries would correspond to the 
broken point group symmetry of the underlying lattice. 
It has been suggested that some of these electronic smec- 
tic and nematic phases may play an important role in the 
physics of strongly correlated electron systems 01 • While 
one may try to understand these phases by following the 
studies of structural liquid crystals, namely introducing 
topological defects in the smectic and nematic phases, it 
turns out that a different approach may be more useful 
in the case of electronic systems. This amounts to study- 
ing the Fermi surface instability of electrons upon tuning 
the interaction strength in a novel channel and/or other 
parameters of the system 0, . 

An effective model in this spirit that incorporates 
the two-body interaction between quadrupolar density of 
electrons has been proposed and analyzed 0, la S B ■ 
It was shown that, in the lattice, the transition to the ne- 
matic phase as a function of chemical potential is gener- 
ically first-order 0| . This leads to abrupt changes in 



various physical properties including the Fermi-surface- 
topology change from a closed to an open Fermi surface 
and Hall constant the other hand, there have 

been several studies of extended Hubbard models that re- 
veal the existence of nematic phases in some parameter 
space of the models |^ ^3 ^3 ■ The direct connection 
between two approaches is not known, but as will be 
mentioned later the effective quadrupolar model seems 
to be consistent with the conditions under which the ne- 
matic phases are seen in the studies of extended Hubbard 
models. 

In this paper, we consider itinerant metamagnetism 
arising from the formation of a nematic phase. Metamag- 
netism in metals refers to a jump in magnetization from 
a low magnetization state to a high magnetization state 
as an external magnetic field increases. The transition 
is first order and associated with a jump in the relative 
volume of the spin-up and spin-down Fermi surfaces. In 
principle, this transition does not have to involve spon- 
taneous breaking of the discrete rotational symmetry of 
the lattice. Here we show that metamagnetic transition 
in metals can occur due to the formation of a nematic 
phase where only the spin-up (or spin-down) Fermi sur- 
face spontaneously breaks the discrete rotational symme- 
try. 

We study a generalized effective Hamiltonian on the 
square lattice with the quadrupolar-density interaction 
including the spin degree of freedom. As the magnetic 
field increases, the Fermi surface volume of the spin-up 
Fermi surface increases and when the spin-up Fermi sur- 
face gets close to the van-Hove singularity, the nematic 
distortion of the spin-up Fermi surface occurs via a first- 
order transition. This leads to the abrupt change in the 
spin-up Fermi-surface topology from a closed to an open 
Fermi surface while the spin-down Fermi surface changes 
continuously. As emphasized in Ref this is an ex- 
ample where the so-called Lifshitz transition of non- 
interacting systems, which involves a continuous quan- 
tum transition between two states with topologically dis- 
tinct Fermi-surface structures without breaking any lat- 
tice symmetry, is avoided in interacting systems. 



2 



When the magnetic field increases further, another 
first-order transition occurs; the spin-up Fermi surface 
changes from an open to a closed Fermi surface. The 
discrete rotational symmetry is now restored in the high- 
field phase. It is shown that two consecutive first-order 
transitions accompany the jumps in the magnetization, 
leading to two metamagnetic transitions. The transition 
turns into second order j3] at some finite temperature 
and the nematic phase is bounded above by second order 
transition. This behavior is reminiscent of that suggested 
for the unusual phase near the metam agne tic quantum 
critical end point of Sr3Ru207 0, 0, 0| and we will 
discuss possible application of our model to this material. 

Effective Hamiltonian; The effective Hamiltonian 
is based on symmetry considerations such that the re- 
sulting order parameter theory has the correct symme- 
tries. In the case of classical liquid crystals, the directions 
parallel and anti-parallel to the rod-like molecule do not 
make difference. As such, the nematic order parameter 
can be represented as a quadrupolar (second order sym- 
metric traceless) tensor built from the spatial directions. 
In two-dimensions, it changes sign under a rotation by 
7r/2 and is invariant under a rotation by tt. In similar 
spirit, one can construct a quadrupolar order parameter 
for electronic systems using the momentum operators of 
electrons, i.e., Qij = piPj — ■^p^Sij. The attractive in- 
teraction between quadrupolar densities will lead to the 
alignment or anti-alignment of electronic momenta (they 
arc equivalent just like the case of classical liquid crys- 
tals), namely the formation of the nematic order. 

The simplest generalization of this idea to the electrons 
with spins on the square lattice leads to the following 
effective Hamiltonian. We also introduce the external 
magnetic field, i?ext, so that the effective Hamiltonian 
only has the Ising degree of freedom in the spin space. 

(ek-Ai)cLcka -^^(c^T^kT -Cki^ki) 

k,Q = Ta k 

+ E F2{q){Tv[Q^{(i)Qahq)]+h.c.} , (1) 

q,a=T,i 

where h — /isiJcxt and /is is the Bohr magneton. Here ek 
is the single particle dispersion and /i is the chemical po- 
tential. On the square lattice, et = —2t{coskx + cos ky). 
As discussed in Ref. , the inclusion of the next-nearest 
hopping does not change the qualitative picture, so we 
will not include this term for simplicity. (5(q) is the 
lattice analog of the quadrupolar density tensor Q con- 
structed from the momentum operators and has the fol- 
lowing form in the square lattice. 

/ cos kx — cos ky sin kx sin ky \ 
k+qa I sinfc^sinfcj^ cos ky — cos k^) ^^"^ 



<3a(q) = Xl"=k 



Here the interaction, i^2(q), can be any arbitrary short- 
ranged interaction such that -F2(q) — > — < in the 
q — > limit. 



Some remarks on this model are in order. In the 
studies of extended Hubbard models, it was found that 
the nematic order (or Pomeranchuk instability) arises 
when there exist attractive (repulsive) interaction be- 
tween electrons in the same (different) patches via for- 
ward scattering 0. This is consistent with our ef- 
fective Hamiltonian where the form factor of the inter- 
action, —F2{coskx — cosky){cosk'x — cosk'y) provides an 
attractive interaction between electrons from the same 
region near (±7r, 0), and a repulsive interaction between 
an electron near (±7r,0) and another near (0, ±7r). Thus 
it appears that our effective Hamiltonian may be related 
to the low energy limit of extended Hubbard models in 
some parameter regime. 

Mean Field Theory; It is sufficient to concentrate 
on the q limit of the interaction by letting F2(q) 
—F2^q,o with F2 > 0. In order to obtain mean field the- 
ory, the quartic interaction is decoupled via the following 
order parameters. 



A„=F2(Q?i(0)), A'„=^^2(Q?2(0)) 



(3) 



The effective single particle dispersion for each spin de- 
gree of freedom now becomes 

e^jj = ~2t{coskx+cosky)+Aa{coskx~cosky)~A'^smkxSm.ky 

(4) 

When a non-zero solution for and A'^ exists, the spin- 
a Fermi surface will spontaneously break the discrete ro- 
tational symmetry of the lattice. It turns out that is 
always zero in our simple model so that the deformation 
of the Fermi surface occurs along the symmetric axes of 
the underlying lattice. Thus we will set = from 
now on. 

The grand canonical free energy density in the mean 
field theory, F{fi, /i, A|, A|) = F| + Fi, is given by 



F„(/i„,A„) = ^^ + FO(/.„,A„) , 
where = ^ + h and /i| = /i — h. Here F^ is 



(5) 



FO(Ai„,A„) = -T J de„P(e„)ln(l+e-(^°-''")/^), (6) 

where V{ea) is the density of states of the renormalized 
single particle dispersion relation for each spin degree of 
freedom and has the following functional form. 



I?(ea) = iVoRe 



1 



.(2) 



-.K 1 



1 - 



it , 



(2t)2 - (^) 



(7) 

where K{x) is the elliptic integral of the first kind and 
A'o = 1/ (27r^t) . Note the persistence of the logarithmic 
van Hove singularities at (eQ,/2)^ = A\ that are remnants 
of the van Hove singularity of non-interacting systems. 
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This leads to non-analytic behavior of the free energy at 
(/.T/2f = A2or(Mi/2)2 = A2. 

As the external magnetic field increases, fj,-\ in- 
creases (decreases) and when the spin-up Fermi sur- 
face gets close to the first Brillouin zone boundary, the 
non-analytic behavior of F-^ near the van Hove filling, 
(^l/2)^ = Ap leads to the first order transition to the 
nematic phase or the abrupt deformation of the spin-up 
Fermi surface. Here the spin-down Fermi surface is far 
away from the van Hove filling so that the spin-down 
Fermi surface only changes gradually. 

Nematic Order and Magnetization; At the zero 
temperature, the free energy, F = F-^ + can be fur- 
ther simplified when A^, ^q,/2 < 2t and -F^ for each spin 
degree of freedom has the following form 




0.1 0.2 h 



FIG. 1: The nematic order parameter, A^, as a function of 
Zeeman magnetic field, h, for F2N0 = 0.1 and jj, = —0.16 in 
unit of 2t at r = 
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-A'a) (8) 



where all energy scales are in units of 2t for simplicity 
and Aq -independent terms are dropped. 

The free energy is investigated function of the 

applied magnetic field, h, for a given chemical potential 
/i. We start with a closed Fermi surface. As the magnetic 
field, h, increases, the spin-up (spin-down) Fermi surface 
volume increases (decreases). For small magnetic fields, 
there is no spontaneous breaking of the lattice symmetry 
in the Fermi surfaces. When h reaches the value that 
makes the spin-up Fermi surface gets close to the van 
Hove singularity, the nematic order parameter for the 
spin-up electrons jumps to a finite value, A| 7^ 0. This 
represents a first order transition from a closed to an open 
spin-up Fermi surface. As h further increases, another 
first order transition occurs such that A| abruptly jumps 
down to zero and the lattice symmetry is restored. To 
summarize, a nematic phase exists in /ii < /i < /12 and is 
bounded by two first order transitions from (to) the low 
(high) field 'isotropic' phases. 

The behavior of the nematic order parameter for the 
spin-up Fermi surface as a function of h is shown in 
Fig.l. Here we choose n/{2t) = -0.16 and F2N0 = 0.1. 
Notice the order parameter jumps to a finite value at 
hi = 0.0428 and jumps down to zero at = 0.277. The 
behavior of the magnetization, M = M-B Sk(('^kT''kT) ~ 
(c^j^Ckx)), as a function of magnetic field is shown in 
Fig. 2. Note that the magnetization jumps up at hi and 
/i2. Thus there exist two consecutive metamagnetic tran- 
sitions. The corresponding Fermi surfaces in three differ- 
ent regions are shown in the inset in Fig. 2. In principle, 
the nematic phase in hi < h < h2 will exhibit anisotropic 
transport properties. In real systems, however, there will 
be domains with two possible orientations of the open 
Fermi surface for spin-up electrons. 

Possible Applications to Experiments; Metam- 




FIG. 2: The magnetization, M, as a function of Zeeman mag- 
netic field for F2N0 = 0.1 and — —0.16. Two sharp metam- 
agnetic transitions are due to the phase transitions between 
the 'isotropic', (a) and (c), phases and the nematic state, (b). 
The inset shows the up-spin and down-spin Fermi surfaces in 
the nematic and 'isotropic' phases. Note that the down-spin 
Fermi surface volume gradually decreases as the Zeeman mag- 
netic field increases, while the up-spin Fermi surface changes 
abruptly at the nematic-'isotropic' transitions. 



agnetism occurs quite often in complex materials, but 
recent experiments on the bilayer ruthenate, Sr3Ru2 07, 
seem to provide rather unique opportunity to look at the 
interplay between quantum fluctuations and metamag- 
netism [lj;ll3- Early experiments on Sr3Ru2 07 suggest 
that the metamagnetic critical end-point can be pushed 
down to zero temperature by changing the direction of 
the magnetic field In this case, one expects that 

the fluctuations near the critical end point is inherently 
quantum mechanical [iH ll^ . Indeed many signatures 
of quantum critical behavior are seen at finite tempera- 
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tures near the putative metamagnetic quantum critical 

end point [11 [ill- 

More recent experiments ,18| on much cleaner samples, 
however, show that the system avoids the zero temper- 
ature quantum critical end point when the angle of the 
magnetic field gets very close to the 'critical' value. In- 
stead of sustaining large quantum fluctuations, the sys- 
tem undergoes two consecutive first order transitions and 
the magnetization jumps at each transition 18] . It has 
been recently suggested that a magnetic field can 
tune an itinerant system so that the majority-spin Fermi 
surface lies close to a van Hove point which may lead to 
a spin-dependent Pomeranchuk instability |2fll | which is 
destroyed by further application of the magnetic field. 

At the phenomenological level, our simple model gives 
rise to a similar behavior. In the experiment, hi and 
/i2 are quite close while our mean field theory gives a 
relatively wide region of the ncmatic phase. However, 
the window of the nematic phase decreases as the value 
of F2N0 decreases. This quantitative feature can also 
change if one considers a different band structure or intro- 
ducing an explicit interaction between spin-up and spin- 
down electrons. On the other hand, we need an attrac- 
tive quadrupolar density interaction to get the nematic 
phase in a range of magnetic fields while the origin of 
this interaction in real system is not clear. Thus it is not 
easy to make a direct connection between our model and 
the experimental results on Sr3Ru207. One can imag- 
ine, however, that when the metamagnetic fluctuations 
become very large, the Fermi surfaces get soft 14] and the 
dominant term in the low energy effective Hamiltonian 
might look like what we have studied. 

Summary and Discussion; In this paper, we con- 
sider a magnetic-field-tuned transition at zero tempera- 
ture between an 'isotropic' (up to the discrete rotational 
symmetry of the lattice) and a nematic phase where 
the spin-up (or spin-down) Fermi surface spontaneously 
breaks the rotational symmetry of the lattice. Such a 
transition on the lattice is generically first order and nat- 
urally gives rise to metamagnetism with the jump in the 
magnetization across the 'isotropic'-ncmatic transition. 
The first order transition turns into second order at some 
finite temperature so that the nematic phase is bounded 
above by second order transition _^]- 

More specifically, we consider a simple model of a 
quadrupolar density interaction between electrons in Zee- 
man magnetic field. As the magnetic field increases, the 
Fermi surface volume of the spin-up (spin-down) elec- 
trons increases (decreases) . When the spin- up Fermi sur- 
face gets close to the van Hove filling, a first order tran- 
sition to a nematic phase occurs. In the nematic phase, 
the spin-up Fermi surface spontaneously breaks the ro- 
tational symmetry of the lattice and becomes an open 
Fermi surface. As the magnetic field increases further. 



another first order transition happens. In the high field 
phase, the lattice symmetry is restored and the spin-up 
Fermi surface becomes closed again. It is shown that 
the magnetization jumps at each transition between ne- 
matic and isotropic phase. Note that, in this case, the 
occurrence of the nematic order is the source of the meta- 
magnetism. 

In order to obtain these results, we assumed an at- 
tractive interaction between the quadrupolar densities of 
electrons. How this interaction can be generated in a 
more microscopic model is an interesting question and a 
subject of future study. 
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